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ABSTRACT 

We present scaling relations between jet power and radio power measured using the Giant Metrewave Radio 
Telescope (GMRT), Chandra and XMM-Newton, for a sample of 9 galaxy groups combined with the Birzan et 
al. sample of clusters. Cavity power is used as a proxy for mechanical jet power. Radio power is measured 
at 235 MHz and 1.4 GHz, and the integrated 10 MHz- 10 GHz radio luminosity is estimated from the GMRT 
610-235 MHz spectral index. The use of consistently analysed, high resolution low-frequency radio data from 
a single observatory makes the radio powers for the groups more reliable than those used by previous studies, 
and the combined sample covers 6-7 decades in radio power and 5 decades in cavity power. We find a relation 
of the form Pj e t oc L~"J for integrated radio luminosity, with a total scatter of CTLrad=0.63 and an intrinsic scatter 
°f cr i,Lrad=0.59. A similar relation is found for 235 MHz power, but a slightly flatter relation with greater scatter 
is found for 1.4 GHz power, suggesting that low-frequency or broad band radio measurements are superior jet 
power indicators. We find our low-frequency relations to be in good agreement with previous observational 
results. Comparison with jet models shows reasonable agreement, which may be improved if radio sources 
have a significant low-energy electron population. We consider possible factors which could bias our results 
or render them more uncertain, and find that correcting for such factors in those groups we are able to study in 
detail leads to a flattening of the Pj e t:L ra di relation. 

Subject headings: galaxies: clusters: general — cooling flows — X-rays: galaxies: clusters — galaxies: active 



1. INTRODUCTION 

X-ray observations of clusters and groups over the past 
decade have provided strong evidence that, despite central 
cooling times signific antly shorter than the Hubble time (e.g., 
ISanderson et al. 12006b. relatively litt l e gas actually cools be- 
low ~0.5 keV dPeterson et al.l l2003t iKaastra et~al]|2004h . It 
is now widely accepted that active galactic nuclei (AGN) in 
the central dominant ellipticals of these systems can reheat 
the gas through a variety of mech anisms dPeterson & F abian 
2006; McNamara & Nulsen 2007, and references therein). 

The most commonly observed evidence of interaction be- 
tween AGN and the intra-cluster or intra-group medium 
(IGM) is the presence of cavities, formed when AGN jets in- 
flate radio lobes, and identified from the resulting X-ray sur- 
face brightness decrement. Cavities provide a relatively sim- 
ple method for estimating the power output of the jets, since 
the mechanical energy required to expel the IGM can be esti- 
mated from the cavity volume and surrounding pressure, and 
the timescale over which the cavity has formed can be esti- 
mated from dynamical arguments. Cavities are expected to 
heat the surrounding gas via the turbule nt wake produced a s 
they rise buoyantly through the IGM (Churazov et al. 2001). 
Studies of cavities have shown that the energies required 
to create them are sufficient to suppre ss cooling in sys- 
tems across a wide range of mass scales dBrrzan et al.1120041 
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Dun netal]l2005l: iRaffertv et al.ll2006l iDunn et al]l20Tol) . 

The relationship between the mechanical power and radio 
emission of AGN jets and lobes is of interest for two main 
reasons: Firstly because it provides in sight into the physi- 
cal nature of the jet (e.g.. IWillott et al.l . 1999): Secondly be- 
cause it allows estimation of the energy available from AGN 
based on more easi l y acq uired radio data (e.g., iBest et all 
120071) . iBirzan et al.l (120041) determined the relation between 
cavity power and 1 .4 GHz radio power, using a sample dom- 
inated by galaxy clusters. However, many cavities are un- 
detected at 1.4 GHz since radiative aging will cause higher 
frequency emission to fade fast est once the j ets c ease to in- 
ject new plasma into the lobes. IBirzan et al.l (120081 hereafter 
B08) addressed the problem of these ghost cavities by mea- 
suring the relation at 327 MHz, and using an estimate of the 
integrated 10 MHz-10 GHz radio luminosity, both of which 
should be more reliable, and both of which produced steeper 
relations. ICavagnolo et al.l (I2010L hereafter C10) extended the 
relation to lower jet and radio powers by adding 21 giant el- 
lipticals to the B08 sample, again finding a steeper slope, but 
were hampered by the poor quality of available archival low- 
frequency radio measurements. 

We have compiled a sample of 18 galaxy groups, chosen to 
show signs of AGN/IGM interactions, and observed both by 
the Giant Metrewave Radio Telescope (GMRT) an d by Chan- 
dra and/or XMM-Newton dGiacintucci et alj|201 ll hereafter 
Gl 1). Of these groups, nine have cavities (identified as decre- 
ments in surface brightness and in some cases as temperature 
or abundance features) which are clearly correlated with radio 
structures, and in this paper we add these to the B08 sample to 
examine the relations between jet mechanical power and ra- 
dio power. Our sample has several advantages over previous 
studies: 1) Our low-frequency data were acquired from a sin- 
gle observatory and are analysed uniformly, making both flux 
density and spectral index measurements more reliable than 
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is possible for data collected from mixed archival sources; 2) 
We are able, for the first time, to measure the integrated radio 
luminosity for a significant number of low-radio-luminosity 
systems, as well as single-frequency powers; 3) We have a 
closer correlation between radio and X-ray morphologies in 
most cases, than is possible either at 1.4 GHz (where many 
cavities are undetected) or with low-resolution low-frequency 
observations (where unrelated sources may be difficult to re- 
move); 4) Our groups have low radio and cavity powers (typi- 
cally P cav < 10 44 erg s _1 and P1400 < 10 24 W Hz" 1 ), filling in a 
region of parameter space sparsely populated in the B08 sam- 
ple. 

We describe our sample and analysis techniques in §|2] Our 
results and their relation to previous work are discussed in 
and we present our conclusions in §|4] A ACDM cosmol- 
ogy with H[)=70 km s" 1 Mpc -1 and £1m=1-^a=0.3 is adopted 
throughout the paper. The radio spectral index a is defined as 
S u oc v~ a , where S„ is the flux density at frequency v. 

2. OBSERVATIONS AND DATA ANALYSIS 

The nine groups in our study were selected to have cavi- 
ties and X-ray data of sufficient quality to allow a reliable 
determination of the properties of the IGM. Table Q] lists the 
properties of the sample. Images and an in-depth description 
of the radio properties of the groups are presented in Gil. 
We do not include the giant sources NGC 315, NGC 383 and 
NGC 7626, since the available data do not allow clear iden- 
tification of cavities. NGC 315 and NGC 383 extend outside 
the XMM field of vi ew, and NGC 7626 is in a merging group 
(Rand aTTet alJ l2009) where determination of IGM properties 
and gravitating mass is unreliable. NGC 140 7 is not in- 
cluded since the cavity tentatively identified by iDong et aT] 
(2010) is on a much smaller scale than the radio emission, 
and its identification and size vary with image processing. 
NGC 741, where a ghost cavity has been previously identified 
dJefha et al.ll2008h is excluded since the visible radio emission 
appears to be unrelated to the cavity, and is contaminated by 
a second AGN within the group. UGC 408 was not listed as a 
cavity system in Gl 1, as no detailed examination of its struc- 
ture has yet been published. However, it is included in this 
sample, since its radio source appears to have cleared an el- 
liptical cavity in the group centre; it was also included in the 
CIO sample. A more detailed description of IGM structures in 
our full sample of 1 8 groups will be presented in a later paper. 

In three cases we had al ready published indiv idual stud- 
ies of the cavities ; AWM 4 dO' Sullivan et alj|2010h. H CG 62 
dGitti etalj|20loh and NGC 5044 dDavid et alj 12009). For 
these groups, we used the gas properties determined from 
these prior analyses. We note that for HCG 62 the gas tem- 
perature and density profiles were extracted from an XMM- 
Newton observation; in all other cases Chandra data were 
used. 

For the remaining systems, we determined the gas prop- 
erties from the longest available Chandra observations. Ta- 
ble [TJ provides basic information on each dataset. Observa- 
tions were downloaded from the Chandra archive and repro- 
cessed using CIAO 4.2 and CALDB 4.3.0 following the stan- 
dar d techniques described in the Chandra analysis thread^ 
and |Q' Sullivan et aT] d20 1 Ot) . Point sources were identified us- 
ing the WAVDETECT task and, with the exception of sources 
coincident with the group-central AGN, removed. As the 
groups typically fill the field of view, background spectra were 

6 http://asc.harvard.edu/ciao/threads/index.html 



drawn from the standard set of ACIS blank sky background 
events files, scaled to produce the same 9.5-12.0 keV count 
rate as that in the target observation. Very faint mode screen- 
ing was applied to target and background observations where 
appropriate. 

Spectra were extracted fro m circular annuli and a depro- 
jected absorbed APEC model (Sm ith et afl2001l) was fitted to 
each set of spectra using XSPEC vl2.6.0k. The absorbing col- 
umn was fixed at the galact ic value determined from the sur- 
vey of Kal berla et alJ (120051) . Redshifts and adopted distances 
for each group are given in Gl 1. Energies below 0.5 keV and 
above 7 keV were ignored during fitting, so as to minimise 
calibration and background uncertainties. 

Jet powers were assessed using the standard approach of 
assuming that the mechanical power of the jet can be approx- 
imated as the energy of the detected ca vities averaged over 
some timescale (e.g. jBirzan et al.l f2.Q04. B08). We define the 
energy of each cavity to be 4pV, the enthalpy of a cavity filled 
with a relativistic plasma. Cavity sizes are defined by match- 
ing an ellipse to their apparent shape, and assuming a line- 
of-sight depth equal to the minor axis. Uncertainties on the 
volume are estimated by assuming a minimum depth of half 
this value, and a maximum equal to the major axis. The un- 
certainty on the volume is, in most cases, the largest contrib- 
utor to the uncertainty on the energy of each cavity and the 
total uncertainty on the cavity power. Outburst timescales are 
typica lly estimated from the buoyant rise time of the cavities 
(th„„ y . IChurazov et al.ll200l|) . but we also estimate the sonic 
timescale (i.e., the travel time between the AGN and farthest 
point of the cavity at the sound speed) and refill timescale of 
the cavities. Our cavity power estimates assume the buoyancy 
timescale, to ensure that they are directly comparable with 
previous studies. The uncertainty on this estimate includes 
the uncertainties on volume, pressure, and tbuoy These val- 
ues are used when fitting the relations between cavity power 
and radio power. The additional uncertainty associated with 
the other timescale estimates is indicated when we plot these 
relations, but is not included in the fitting process. 

The analysis of our GMRT radio data is described in Gl 1. 
Most of the groups in the sample were observed at both 
235 MHz and 610 MHz, but of the cavity systems included 
in this study, NGC 5813 was only observed at 235 MHz and 
NGC 5846 only at 610 MHz. Radio powers at each frequency, 
P^, are estimated as P ; ,=47rD£(l+z) Q ~ 1 S z ,, where S„ is the 
flux density at frequency v, a the spectral index, z is the red- 
shift and Dl is the luminosity distance to the source. Radio 
powers at 1.4 GHz were also estimated from flux densities 
mainly derived fro m the NRAO VLA Sky Survey (NVSS, 
iCondon et al.l[T998T) . 

Since it is desirable to compare our 235 MHz GMRT mea- 
surements with the 327 MHz VLA radio powers used by B08, 
we have used the spectral indices given in B08 and Gil to 
correct each set of flux densities to the other band, and de- 
termined the best-fitting relationship between cavity power 
and radio power for both frequencies. It should be noted 
that while Gl 1 determined 235-610 MHz spectral indices for 
the majority of their sample, only a lower limit on the spec- 
tral index could be determined for NGC 5044 ( a > 1.9), and 
no index could be determined for NGC 5813 or NGC 5846, 
since they were only observed at one frequency. When cor- 
recting those data to 327 MHz, we adopt a spectral index of 
a=l.9 for NGC 5044, the 235 MHz- 1.4 GHz spectral index, 
a=0.94, for NGC 5813, and the 610 MHz-1.4 GHz spectral 
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Table 1 

Summary of observations and data for each group 



Group 


Instrument 


ObsID 


*exp 


D 

(wipe; 


Kin 

235 


Pl400 

( i n 24 w M7~i *t 

i iu vv n£ f 


P235 
Ifl 24 W T47 _li t 

IO VV Il£ 1 


Lradio 
HO 42 pt-ct e^l"i 


P " 

r cav 
MO 42 pt-ct c - !^ 

dig t> ) 


AWM4 


ACIS-S 


9423 


74.5 


128.7 


0.67 


1.41±0.07 


6.31±0.50 


0.14±0.02 


53 67 +iJj5 


HCG 62 
UGC 408 


ACIS-S 
XMM 
ACIS-S 


10462 
050478501 
11389 


C*H 1 

o/.l 
122.5 
93.9 


57.8 
61.9 


1.15 
0.53 


U.UU21±U.UUU1 

0.83±0.04 


U.UloztU.UUl 

2.57±0.21 


U.UUUzy±U.UUUU3 
0.071 ±0.008 


3 - 49 -1.66 
3.53 


NGC 507 


ACIS-I 


2882 


43.5 


69.3 


1.14 


0.060±0.003 


0.68±0.05 


0.011±0.001 


20.68^^2 


NGC 4636 


ACIS-I 


3926 


74.6 


13.2 


0.61 


0.0017±0.0001 


0.0053±0.0004 


0.00013±0.00002 




NGC 5044 


ACIS-S 


9399 


82.7 


38.2 


>1.90' 


0.0065±0.0003 


0.042±0.003 


0.0019±0.0002 


9 00+2.55 


NGC 5813 


ACIS-S 


9517 


98.8 


26.8 


0.94* 


0.0014±0.0001 


0.0079±0.0006 


0.00013±0.00002 


1 64+0 29 

1 - w -0.20 


NGC 5846 


ACIS-I 


7923 


90.0 


24.3 


0.65 c 


0.0015±0.0001 




0.000 11 ±0.00001 


1 71 +0.50 
'-0.44 


NGC 6269 


ACIS-I 


4972 


39.6 


142.9 


0.73 


0.141 ±0.007 


0.65±0.05 


0.013±0.002 


10 02+ 2 20 


Note. 


" Only a lower limit on the 235-610 MHz spectral index, 


610 

a 235' 


could be estimated for NGC 5044, as extended structures at 235 MHz are undetected 



at higher frequencies (see Gl 1). * NGC 5813 spectral index measured between 235 MHz and 1.4 GHz. c NGC 5846 spectral index measured between 610 MHz 
and 1.4 GHz. d Cavity power estimated based on the buoyancy timescale and including uncertainties on volume, pressure, and timescale. 



index, a=0.65, for NGC 5846. 

We also estimate the integrated radio luminosity of our 
groups in the 10 MHz- 10 GHz band, assuming a powerlaw 
spectrum and using the same spectral indices. These are not 
perfectly comparable with the integrated radio power esti- 
mates of Birzan et al., which are based on model fits to flux 
densities measured at 3-8 frequencies for each system. How- 
ever, the inaccuracy introduced by our use of a simple spec- 
tral index is unlikely to be significant compared to the large 
uncertainties in cavity power. Low-frequency spectral index 
measurements such as ours should not be strongly affected 
by radiative losses unless the break frequency is below our 
lowest measurement frequency (235 MHz). It seems unlikely 
that any of our systems are so old. Two groups in the sam- 
ple, NGC 50 7 and AWM 4, have unusually high radia tive 
age estimates (Murg ia et al.1201 ltlGiacintucci et a l. 2008) but 
break frequencies of 300-450 MHz, above our lower fre- 
quency bound. If the spectral indices are accurate, then we 
are unlikely to over-estimate the steepness of the spectrum 
and the power at low frequencies, as would probably be the 
case if high-frequency spectral indices were used. Detailed 
modelling of the radio spectra for our group sample and re- 
vised estimates of integrated radio power will be presented in 
a forthcoming paper. 

3. RESULTS AND DISCUSSION 

Figure [T] shows cavity power vs. radio power at 1.4 GHz 
and 235 MHz for our nine groups and the 24 systems de- 
scribed by B08. At both frequencies, the scatter among the 
points is considerable. The effect of radiative aging on the 
radio sources causes preferential fading of emission at higher 
frequencies. We therefore expect the low-frequency measure- 
ments to provide a more accurate estimate of the true radio 
power in older systems. Figure [2] shows the relationship be- 
tween cavity power and the integrated radio luminosity in the 
10 MHz- 10 GHz band. The integrated luminosity should be 
a superior measure of the radiative power compared to esti- 
mates at a single frequency, since it accounts for variations in 
spectral index between sources. However, we still see a large 
degree of scatter amongst the data points. 

We used the bivari ate correlated error and int rinsic scat- 
ter (BCES) algorithm (lAkritas & Bershadvl[l996l) to perform 
linear regression fits to the data, determining the best fitting 



power-law relationship between cavity power, P cav , and the 
radio powers at each frequency (P1400, P235, P327) or the in- 
tegrated radio luminosity L ra di . Using the orthogonal BCES 
regression to the parameters in log space, the best-fit relations 
are: 

log P cav = 0.7 1 (±0. 1 1) log P 235 + 1 .26 (±0.12) (1) 
log P cav = 0.72 (±0. 1 1) log P 327 + 1-38 (±0.12) (2) 
log P cav = 0.63 (±0.10) log P 14 oo + 1 -76 (±0.15) (3) 
log P cav = 0.71 (±0.1 1) log L radio + 2.54 (±0.21) (4) 

where P cav and L ra di are in units of 10 42 erg s" 1 , and P235, 
P 327 and P1400 in units of 10 24 W Hz" 1 . 

We estima ted the total scatter in the data using the technique 
described in iPratt et all (J2009), which is based on weighted 
estimates of the orthogonal distance between the data points 
and best fit lines. We find the total scatter for the radio power 
relations to be eri4oo=0.68 dex and (7235=0.62 dex, while the in- 
tegrated radio luminosity relation has a scatter of OLi- a d=0.63. 
The intrinsic scatter in the data can be estimated by subtract- 
ing the contribution expected from the statistical errors. We 
found that the scatter about the radio power relations was 
c;,i400=0.65 dex and (7, 235=0. 58 dex, and the scatter about the 
integrated radio luminosity relation was cr,,Lrad=0.59 dex. 

These relations can be compared with those derived by B08, 
including th e corre ction to the 327 MHz relation given by 
iBrrzan et al.l (120101) and using the integrated radio luminos- 
ity from their sources as a whole rather than from the lobes 
alone: 

log P cav = 0.51 (±0.07) log P 327 + 1 .5 1 (±0. 12) (5) 

log P cav = 0.35 (±0.07) log P1400 + 1 .85 (±0.10) (6) 

log P cav = 0.49 (±0.07) log L radio + 2.32 (±0.09) (7) 

However, these relations were estimated using an ordinary 
least-squares regression, minimising the distance between fit 
line and data points in the P cav axis only. Refitting these re- 
lations using the orthogonal BCES regression steepens them 
and increases the uncertainties, giving slopes of 0.67 ±0.1 9 
at 327 MHz, 0.57±0.17 at 1.4 GHz and 0.68±0.19 for the 
integrated radio power. HCG 62, the only group in the B08 
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Figure 1. Cavity power vs. radio power at 1.4 GHz (left) and 235 MHz (right). Systems in the sample of B08 are marked by grey triangles, members of our 
groups sample by black circles, la uncertainties on radio power and cavity power (calculated using the buoyant rise time of the cavities) are indicated by error 
bars. In many cases the uncertainty on the radio luminosity is smaller than the data point. For our groups, additional narrow-width error bars indicate the 1 cr 
uncertainty range allowing for alternate measures of cavity age (the sonic and refill timescales). The Birzan et al. point with the lowest radio power is HCG 62, 
which is also included in our GMRT groups sample. This point is included for comparison, but was excluded from our analysis. The solid fit line indicates our 
BCES regression fit to the data points. The dotted line indicates the relation found by B08, the dashed line the relation found by CIO. For the 235 MHz relation, 
we have used the spectral indices given by B08 to correct their 327 MHz data to 235 MHz, so as to allow a direct comparison with our data. The normalizations 
of the Birzan or Cavagnolo fit lines are not corrected for frequency differences, but the lines are indicative of the relative gradient of the different fits. 




10" 4 10" 3 10" 2 10-' 10° 10 1 10 2 10 3 



lo g( L „<uc, [ 1042 er 8 s_1 ]) 

Figure 2. Cavity power vs. integrated 10 MHz- 10 GHz radio power. Sym- 
bols are the same as in Fig[T] 

sample, has a strong influence. Excluding this system, fits 
to the B08 dataset have gradients of 0.81±0.30 at 327 MHz, 
0.62±0.28 at 1.4 GHz, and 0.78±0.30 for the integrated ra- 
dio power. This steepening, and the large uncertainties on the 
gradients, emphasize the need for a large dynamic range if 
gradients are to be measured accurately from data with such a 
large degree of scatter. 



B08 also found steeper relations when fitting only those 
systems where the radio emission fills the cavities seen in the 
X-ray, excluding systems with ghost cavities. A comparison 
of radio and cavity power on this basis is likely to be more re- 
liable, since the radio emission is clearly associated with the 
lobes which inflated the cavity. Our relations have gradients 
which are steeper than the published B08 relations (though 
consistent with them at the 1.6-2.3ct significance level), and 
in agreement with those derived from the Birzan et al. sam- 
ple using orthogonal regression. The total scatter about the 
B08 radio luminosity relation is ~0.83 dex, suggesting that 
o ur groups are mor e close ly clustered about the relation. 

Cavagnolo et al. ( 2010), using orthogonal BCES regres- 
sion, found the following relations between radio power and 
cavity power: 

log P cav = 0.64 (±0.09) log P200-400+ 1.54 (±0.12) (8) 

log P cav = 0.75 (±0.14) log P 14 oo + 1 -91 (±0.18) (9) 

Unsurprisingly, given the large overlap between our sam- 
ple and that of C10 (The C10 sample of 21 ellipticals in- 
cludes 8 of our 9 groups, and both sets of fits have the B08 
sample in common) their best-fitting relations are very simi- 
lar to those we derive, with gradients in agreement to within 
the lcr uncertainties at both 1.4 GHz and at lower frequen- 
cies. C10 drew their low frequency radio mea surements from 
the CATS database ( Verkhoda nov et al.ll 19971) . a compilation 
of data from numerous radio catalogues. We expect there to 
be significant variations in sensitivity and resolution among 
measurements from different observatories and surveys. Our 
GMRT 235 MHz were observed and processed uniformly, and 
we therefore expect our data points to be more reliable. Com- 
parison of our 235 MHz radio powers with the 200-400 MHz 
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powers of CIO suggests that agreement is fairly good for the 
most radio luminous systems (P235 > 10 23 W Hz -1 ). How- 
ever, for the fainter sources, disagreements of a factor of ~10 
are typical. This may be explained by issues such as the inclu- 
sion of unresolved background sources and the non-detection 
of faint structures. The total scatter about the CIO relations is 
comparable with ours at low frequencies (cr2oo-400=0.61) and 
slightly larger at 1.4 GHz (cri 4 oo=0.78). 

3.1. Comparison with theoretical models 

IWillott et al.l |l999) determined a theoretical relation be- 
tween radio luminosity and jet power which is now widely 
used to estimate the kinetic energy output o f AGN (e.g., 
lHardcastle et al.l 120071: ICattaneo & Bestl 120091) . This model 
takes the form 

P jet = 3 x 10 38 / 3/2 P^ W, (10) 

where Pj et is the total jet power, P151 is its observed 

151 MHz radio power in units of 10 28 W Hz" 1 Sr~\ and the 
factor / includes a variety of unknown factors affecting the 
normalisation of the relation. Under the assumption that the 
spectra of radio sources are simple powerlaws, the normali- 
sation of the model is dependent on frequency, but its gradi- 
ent is not. Since P cav is a proxy for the total jet power, it is 
straightforward to compare the model gradient with those of 
measured Pcav^Pradio relations. C10 compared their best fitting 
Pcav:Pi400 relation with the model, finding that its gradient of 
^0.86 agreed with their measured gradient within errors. At 
lower frequencies, which should be more reliable, the model 
agrees with the C10 P C av:P200-400 relation within ~2.5cx con- 
fidence bounds. The normalisation of the model is dependent 
on a number of factors, including k, the ratio of energy in non- 
radiating particles to relativistic electrons. C10 found that k 
values of tens to thousands were required to bring the model 
normalisation into agreement with their data. Ignoring their 
relative normalisations, the gradient of our 235 MHz relation 
would agree with the model within 1.5cr, and our 1.4 GHz 
relatio n within 2.5a. 

The Wil lott et al.l d 19991) model is based on relatively 
straightforward synchrotron physics, but includes a number 
of important assumptions. The jet power is proportional 
to the energy density of the relativistic plasma assuming a 
minimum energy magnetic field, u me , following the relation 
3/2 3 

Pj et oc u„ie oc B me , where B me is the minimum energy mag- 
netic field. Willott et al. define B me in terms of observabl e 
radio properties following the method given bv lMlievl (fl980). 

which can be simplified to B me oq(v^ 2 a — v X J 2 a ) 2 / 1 7 where 
v\ and vi are frequencies defining the observable radio band 
(10 MHz- 100 GHz for Willott et al.) and a is the spectral in- 
dex. However, this assumes that the electron population in the 
radio-emitting plasma only contains particles with energies 
such that they radiate in this observable band, and Willott et 
al. note that the jet power will depend critically on the choice 
of low frequency cutoff. 

Electrons with energies too low to produce observable 
emission may contribute a significant fraction of the energy in 
radio jets and lobes, since if they are present they are likely to 
make up the majority of the electron population. Many obser- 
vational studies define the range of electron energies in terms 
of a range of Lorentz factors, typically with a low energy cut- 
off of 7 m in=10 or 100. In a magnetic field of ^1/iG, such elec- 
trons will produce synchrotron radiation at ^0.5-50 kHz, well 



below the observable band. IWorrall & Birkinshawl (120061) 
provide a relation between B me and observable radio parame- 
ters using Lorentz factors rather than frequencies, which takes 
the formB„, e oc (7™ 2a -7i7 2a ) 1/(a+3) . If we use this definition 
to determine u me and therefore Pj et , the dependence of the gra- 
dient of the Willott et al. model on the spectral index becomes 
clear, giving a relation of the form: 



This reduces to a gradient of 6/7 for a = 0.5. However, the 
measured 610-235 MHz spectral indices of the Gil sample 
cover the range 0.53-1.44. Excluding systems for which only 
limits are available, the mean spectral index is o$%$ = 0.95, 
implying a Pcav^Pradio gradient of 0.76. This is considerably 
closer to our best measured gradients of ~0.7 for 235 MHz 
power and integrated radio luminosity. A spectral index of 
0.8, often used as a typical value for extragalactic sources 
(ICondonlll992h would give a model gradient of 0.79. Our 
best fitting P C av:L ra di relation gradient of 0.71 would suggest 
a~1.2. 

Willott et al. developed their model to examine the rela- 
tionship between low-frequency radio power and narrow line 
optical luminosity in a sample of 7C and 3CRR radio sources. 
They found no correlation between 151 MHz spectral index 
and the residual from their best fitting radio to optical rela- 
tion. Such a correlation would be expected if jet power varied 
strongly with spectral index. They argue that this suggests 
that low-energy electrons make only a minimal contribution 
to jet power. However, the scatter in their dataset is large (as 
it is in ours) and their sample contains a wide variety of source 
types, for which spectral index may be dominated by emission 
from different physical regions. FR-II radio galaxies are also 
common in their sample, while our groups host only moder- 
ately powerful FR-I sources and a few FR-I/FR-II transition 
systems. It therefore seems possible that the slope of the re- 
lation could be steeper. Such a change in gradient may have 
implications for models which have used the Willott et al. re- 
lation to estimate the energy output of the population of AGN 
from their radio luminosity function, with a shallower gradi- 
ent implying greater mechanical power available from lower 
luminosity jets. 

3.2. Uncertainties and potential biases 

Several factors could affect our estimates of the 4pV cavity 
power, or the accuracy of these estimates as a proxy for the 
mechanical power of the jets. These include: 

• Cavity volume. Identification and characterisation of 
the cavity is inherently subjective, and dependent on the 
quality of the X-ray data, the angular and physical sizes 
of the cavity, its position in the group and other factors. 
Using HCG 62, which has several Chandra and XMM- 
Newton datasets with a range of exposures available, 
we tested the effects of performing independent spectral 
deprojections and having different researchers estimate 
cavity size. We found differences between estimates 
of up to a factor of 2 (~0.3 dex) in total cavity power, 
despite the relatively simple morphology of the cavities. 

• Cavities at large radii are significantly more difficult to 
identify, owing to the decrease in X-ray intensity with 
radius. In both HCG 62 and NGC 5044, 235 MHz ra- 
dio maps reveal lobes at large radii, beyond the cavities 
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identi fied in the X-ray dDavid et all 120091: iGitti et all 
l2Ol0h . In both cases there is some evidence of the pres- 
ence of a cavity in the X-ray (e.g., a low abundance 
region suggesting multi-phas e gas coincident w ith the 
detached lobe in NGC 5044 iDavid et al.ll2011l) . and if 
such cavities exist and have sizes similar to the radio 
lobes, this would increase our estimate of the AGN 
power by a factor of ^2-10. Multiple cavity pairs in 
individual groups also give some idea of the range of 
cavity powers. In both NGC 5044 and NGC 5813, cav- 
ity powers for individual cavities vary by factors of up 
to 10. 

Very old or young cavities are unlikely to be detected. 
Young, small cavities would have a minimal X-ray sur- 
face brightness decrement. Their apparent cavity pow- 
ers would also be small, but since they would likely be 
highly over-pressured, cavity power would underesti- 
mate jet power. Very old cavities would be expected 
to have risen to very large radii or to have already bro- 
ken up into smaller structures; in either case the surface 
brightness decrements would be small. Their cavity 
powers would be low, owing to their long timescales. 
Non-detection of such cavities in groups could bias the 
Pcav^Pradio relations to steeper gradients, but it is unclear 
how common such cavities are, or whether similar bi- 
ases might also affect galaxy clusters. 

Shocks driven by the AGN outburst may contain a large 
fraction of the energy released. Inclusi on of the shocks 
detected in HCG 6 2 and NGC 5813 (IGitti et alJl2Q10b 
iRandall et alJl201 H) would increase our estimate of the 
AGN power output by a factor ^10. Given the diffi- 
culty of detecting shocks (deep Chandra observations 
are generally required) we cannot know whether we 
are missing a significant energy contribution in other 
groups. 

Uncertainties in outburst timescale will also affect 
power estimates. We have used the buoyancy timescale 
and its uncertainties when fitting the P cav :Pradio rela- 
tions, and included additional error bars to show the 
larger uncertainty range associated with different dy- 
namical timescale estimates, but even these may not 
be accurate, particularly for systems which have ex- 
panded supersonic ally for a significant fraction of their 
lifespan (e.g., UGC 408). Very old sources may also 
provide poor estimates if the buoyant velocity over- 
estimates their true rate of rise; radiative ages based 
on synchrotron losses can exceed the the dynami- 
cal timesc ales in such systems (e.g., in AWM 4 an d 
NGC 507. IP' Sullivan et alJMTot iMurgia et al.ll20ll . 

Filling factors of less than unity for radio lobes could 
also render volume estimates for cavities inaccurate, af- 
fecting both the energy and dynamical timescales of 
the outburst. Our study of AWM 4 suggested that 
the lobes may have a filling factor as low as 0=0.2 
dO' Sullivan et aUl20T(ih . 

Jet orientation could affect our estimate of the position 
of cavities in the IGM, with jets close to the line of sight 
producing cavities which appear to be at smaller radii 
than is the case. This will lead to underestimates of cav- 
ity enthalpy and outburst timescale. Simulations sug- 



gest that cavity powe rs will typically be with in a factor 
of 3 of the true value (Me ndygral et al.ll20TTl) . 

• AGN-driven "weather", turbulent motions in the IGM 
induced by AGN jets, could affect the position of cav- 
ities. NGC 5044 may provide an example, with sev- 
eral s mall cavities foun d at similar radii in the group 
core (Dav id et al.11200 90. Timescales for such cavities 
are probably underestimated, leading to overestimates 
of cavity power. 

For our sample of groups, it appears that the effect of such 
biases may be to steepen the relation between cavity power 
and radio power. Correcting the cavity power estimates for 
HCG 62, NGC 5044, NGC 5813 and AWM 4 to include con- 
tributions from the probable outer cavities and shocks, and the 
low filling factor of the AWM 4 radio lobes, we find that the 
gradient of the P cav :L rac iio relation flattens to 0.62±0.12. This 
is still consistent, within uncertainties, with our initial fit, but 
we note that these are perhaps the four most carefully studied 
systems in our sample, and deeper observations of other sys- 
tems would likely lead to similar corrections. For the systems 
we can study in detail, corrections which reduce cavity power 
tend to be small, whereas the effect of including shocks and 
possible additional cavities can increase the estimated AGN 
power by much larger factors. 

The impact of such biases on galaxy clusters is less clear. 
These are generally more X-ray luminous systems with larger 
cores, hosting more powerful radio sources, which might sug- 
gest that cavities will be more easily detected. The intra- 
cluster medium is also likely to be able to confine more pow- 
erful AGN outbursts, which in groups might simply tunnel 
out to large radii where any cavities would be difficult to de- 
tect (e.g., the 'poorly-confined" class of objects in the C10 
sample). However, weak shocks may be more easily detected 
in galaxy groups, since the emission lines produced by gas at 
~1 keV significantly improve our ability to accurately mea- 
sure temperature. Resolution may also be a problem for more 
distant clusters, producing a bias toward the detection of large 
cavities associated with more powerful AGN outbursts. A sta- 
tistical approach to the effects of such biases, based on simu- 
lations, would perhaps provide insight into this problem. 

4. CONCLUSIONS 

We have estimated the relations between mechanical jet 
power and radio luminosity, adding nine groups selected from 
the sample of Gl 1 and observed with the GMRT and Chandra 
or XMM-Newton to the B08 sample, which consists primarily 
of galaxy clusters. We find P cav :Pradio relations with gradi- 
ents of ^0.7 for both the low-frequency radio power (235 or 
327 MHz) and the integrated radio luminosity, with total scat- 
ters of ov,Lrad=0.59 and ct (j 235=0.58 dex. The 1.4 GHz relation 
is somewhat flatter (gradient ^0.6) and has a slightly larger 
scatter, a, ;.i4oo=0.65 dex. In agreement with previous studies, 
this suggests that low-frequency and broad band radio mea- 
surements are superior indicators of cavity power. The weaker 
correlation between 1.4 GHz radio structure and the cavities 
identified from the X-ray makes high frequency emission a 
poor choice for such studies. Our fitted slopes are signifi- 
cantly steeper than those found by B08, but this is unsurpris- 
ing since the B08 relations were determined using a different 
regression technique, which will tend to produce a shallower 
gradient. Using the same BCES orthogonal regression used 
for the fits to the combined dataset brings the B08 relations 
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into agreement with our results. Our P ca v:P235 relation has a 
similar gradient to the P ca v:P200-400 relation of CIO. This is 
expected, given the overlap between samples, but direct com- 
parison of the radio powers suggests that our GMRT measure- 
ments are more reliable for low -power radio galaxies. 

O ur PmviP?.'^ relation is somewhat flatter than the widely 
used IWillott et alj (Q999) model of jet mechanical and radio 
power, though they agree within 1.5a uncertainties. We find 
that inclusion of electrons with low Lorentz factors (which 
cannot be directly observed) could change the gradient of the 
Willott model, making it dependent on the radio spectral in- 
dex. In this case, using the mean spectral index of the Gil 
sample, the model would agree more closely with our ob- 
served relation. A variety of factors could bias or increase 
the uncertainty of our measurements, and we conclude that at 
least for galaxy groups these may have a serious impact on our 
cavity power estimates. Correcting for these factors, in those 
groups where the quality of radio and X-ray data is sufficient 
to allow detailed study, produces a flatter P cav :P, a( jio relation. 
However, it is unclear whether this would be the case in all 
groups, or in galaxy clusters, and simulations of AGN feed- 
back across a wide range of mass scales and outburst powers 
are probably required to resolve this question. 
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